AD-A155  127  PARAMETER  ESTIMATION  IN  COMMUNICATION  SVSTEM  TRACKING  1/1 
SATELLITE  OBSERVATIONS (U)  NAVAL  POSTGRADUATE  SCHOOL 
MONTEREV  CA  V  A  TSAFARAS  DEC  84 


UNCLASSIFIED 


F/G  17/2 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  1 963-A 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


iiilr 


DTIC 

ELECTE 
JUN2  0  1985 


THESIS 


Approved  for  public  release;  distribution  unliraitec 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


SECURITY  CLASSIFICATION  of  this  pace  fWhen  Dele  Entered) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

WSmBmSl 

4.  TITLE  tend  Subtitle) 

Parameter  Estimation  in  Communication 
System  Tracking  Satellite  Observations 

5.  TYPE  OF  REP0Pt  4  PERIOD  COVERED 

Master's  Thesis; 

December  19S4 

5.  PERFORMING  ORG.  REPORT  NUMBER 

7.  auThorc*; 

Vassilios  Ath.  Tsafaras 

8.  CONTRACT  OR  GRANT  SuMflERfiJ 

S.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Postgraduate  School 

Monterey,  California  93943 

10.  program  element,  project,  task 
AREA  8  WORK  *n|T  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  93943 

12.  REPORT  Date 

December  1384 

13.  NUMBER  Cc  3  *  G  SS 

82 

14.  MONITORING  AGENCY  name  4  AOORESSftf  different  Item  Controlling  Office) 

15.  SECURITY  ClASS  ,  ol  true  report) 

UNCLASSIFIED 

15«.  DECLASSIFICATION-  DOWNGRADING  1 

SCHEDULE 

1«.  DISTRIBUTION  STATEMENT  (of  thia  Report) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  the  obttrocl  entered  in  Bloc*  30.  It  different  from  Report) 

18.  supplementary  notes 

If.  KEY  WOROS  r  Continue  on  rmverae  aid*  if  necmaaary  and  Identify  by  block  number) 

Parameter  Estimation;^  Kalman  Filter,  _/ 

/ 

20  ABSTRACT  fConttnua  on  rmvaraa  aide  if  nacaaaary  and  idantity  by  block  number) 

The  estimation  of  parameters  from  a  satellite  communication 
system  is  often  through  the  use  of  Kalman  filtering.  In  this  work 
the  location  of  the  eye  of  a  hurricane  is  estimated  from  satellite 
observations.  A  comparison  with  a  posteriori  meteorologist's 
analysis  was  attempted.  An  adaptive  gating  scheme  was  employed 
in  the  filter  to  accommodate  "manueuvers"  of  the  storm. 

The  observations  were  at  random  intervals  and  also  came  from 
several  different  sources  (aircraft  and  land  based  radar  as  well 

00  i  jan *7 j  1^73  eoition  of  i  nov  «s  is  obsolete  as  satellite)  . 

1  N  0102-  LF-  01 4-  6601  ^  - 


SECURITY  CLASSIFICATION  of  this  PACE  (When  Data  Entered) 


Approved  for  public  release;  distribution  unlimited 


Author: 

Approved 


Parameter  Estimation  in  Communication 
System  Tracking  Satellite  Observations 

by 


Vassilios  Ath.  Tsafaras 
Lieutenant,  Hellenic  Navy 
B.S.,  Greek  Naval  Academy,  1973 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1984 


Vassilios  Ath . _jTsafaras 


by : 


H.A.  Titus,  Thesis  ; 


Advisor 


A.  Gerba,  Jr.,  Second  Reader 


H.  Rigas,  Chairman,  Department  of 
Electrical  and  Computer  Engineering 


J .  Dye 


of  Science  and  Engineering 


ABSTRACT 


The  estimation  of  parameters  from  a  satellite 
communication  system  is  often  through  the  use  of  Kalman 
filtering.  In  this  work  the  location  of  the  eye  of  a 
hurricane  is  estimated  from  satellite  observations.  A 
comparison  with  a  posteriori  meteorologist's  analysis 
was  attempted.  An  adaptive  gating  scheme  was  employed 
in  the  filter  to  accommodate  "manueuvers"  of  the  storm. 

The  observations  were  at  random  intervals  and  also 
came  from  several  different  sources  (aircraft  and  land 
based  radar  as  well  as  satellite) . 


.'...'cession  For 

I.TIS  T':A&I 
P'CIC  TV;1 


TABLE  OF  CONTENTS 

I.  INTRODUCTION - 

II.  KALMAN  FILTERING  TECHNIQUES - 

III.  ERROR  ELLIPSOIDS - 

IV.  SATELLITE  TRACKING-SIMULATION  RESULTS 

V.  RANDOM  TRACKING - 

VI.  CONCLUSIONS - 

APPENDIX  A  -  COMPUTER  ALGORITHM - 

BIBLIOGRAPHY - 

INITIAL  DISTRIBUTION  LIST - 


4 


LIST  OF  TABLES 


1.  OBSERVED  DATA  (SATELLITE  FIXES) -  17 

2.  BEST  TRACK  DATA -  IB 

3.  ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE -  22 

4.  ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 

AFTER  CORRECTION -  28 

5.  BEST  TRACK  DATA  (FICTITIOUS  STORM) -  33 

6.  OBSERVED  DATA  (FICTITIOUS  STORM) -  34 

7.  ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 

(FICTITIOUS  STORM) -  37 


5 


LIST  OF  FIGURES 


1.  Trajectory  of  Nelson  Storm -  19 

2.  Gains  of  Position  and  Velocity -  20 

3.  K.F  Track  and  Prediction  Ahead  in  Latitude -  21 

4.  Latitude  Errors -  23 

5.  Longitude  Errors -  2  4 

6.  Trajectory  of  Nelson-Storm  (with  corrected  algorithm) 25 

7.  Gains  of  Position  and  Velocity  (with  corrected 

algorithm) -  26 

8.  K.F  Track  and  Prediction  Ahead  in  Latitude  (with 

corrected  algorithm) -  27 

9.  Latitude  errors  (with  corrected  algorithm) -  29 

10.  Longitude  Errors  (with  corrected  algorithm) -  30 

11.  Trajectory  of  Fictitious  Storm -  35 

12.  Gains  of  Position  and  Velocity  (fict.  storm) -  36 

13.  Latitude  Errors  (fict.  storm) -  38 

14.  Longitude  Errors  (fict.  storm) -  39 


I .  INTRODUCTION 

Satellite  communications,  especially  digital,  are  well 
under  way.  One  of  the  many  observations  that  satellites 
provide  are  the  meteorological.  These  give  the  location  of 
a  storm  or  a  typhoon  in  terms  of  geological  coordinates.  It 
is  possible  to  have  images  of  the  store  or  typhoon  cloud 
cover . 

These  observations  occur  quite  randomly  in  time. 
Meteorologists  like  to  present  their  forecasts  equally 
spaced  in  time. 

The  attempt,  here,  is  to  try  to  compare  an  adaptive 
Kalman  filtering  algorithm  to  estimate  the  location  of  the 
eye  of  a  storm  with  optimum  track  values  that  are  given  from 
the  meteorologist ' s  analysis.  The  filtering  process  combines 
all  the  available  measurement  data  with  prior  knowledge  about 
the  system.  It  produces  an  estimate  of  the  location  of  the 
storm  (latitude,  longitude)  in  such  a  manner  that  the  mean 
square  error  is  minimised.  The  parameters,  which  are  essential 
design  elements  of  a  Kalman  filter,  are  the  measurement 
noise  covariance  matrix,  R,  the  excitation  covariance,  0, 
the  initialization  covariance  of  error  in  the  filter  itself, 

P  ( 1  / 0 ) ,  and  the  transition  matrix,  0. 

In  most  physical  processes  that  one  desired  to  track, 
many  of  these  parameters  change  during  the  process  of 


tracking.  The  measurement  noise  associated  with  the 
observations  can  change  if  a  different  sensor  is  used,  or 
a  similar  sensor  obtaining  measurements  from  some  different 
geometry  relative  to  the  object  being  tracked.  If  the  object 
being  observed  is  acted  upon  by  external  forces,  then  the  Q 
matrix  should  be  changed  to  account  for  these  external 
excitations.  Most  processes  being  tracked  will  change  in  their 
dynamic  characteristics  during  the  observation  time  and  so  the 
transition  matrix  ideally  should  also  be  changed.  Further, 
the  time  between  observations  quite  often  occurs  randomly 
in  time.  All  of  these  things  bring  about  the  need  to  change 
the  Kalman  filter  to  adapt  as  the  process  changes. 

It  is  possible  to  change  the  parameters  of  the  Kalman 
filter  by  sensing  the  error  between  the  observation  and  the 
prediction  from  the  track  of  what  that  observation  should  be. 

If  a  gate  is  established  representing  95%  of  the  normal 
random  perturbations  of  the  process,  then  when  this  error 
exceeds  the  magnitude  of  the  gate,  one  can  reasonably 
ascertain  that  the  filter  is  no  longer  properly  representing 
the  observed  process.  In  the  work  attacked  here,  real  data 
was  obtained  from  satellite  observations  and  the  qualitative 
observation  error  was  established  (PCN  #) . 

The  error  covariance  matrix  in  terms  of  error  ellipsoids 
along  the  track  gives  a  measurement  for  the  worthiness  of 
the  algorithm  parameters. 
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KALMAN  FILTERING  TECHNIQUES 


In  a  linear,  discrete  system,  the  state  and  measurement 
equations  are  given  by 

x  (k+1)  =  0(k)x(k)  +  £w(k) 

and 

zjk)  =  H(k)x(k)  +  v(k) 

where  x  is  the  state;  _P  is  the  transition  matrix;  _T  is 
excitation  noise  matrix;  H  is  the  measurement  matrix;  w 
and  v  are  the  excitation  and  measurement  noise  correspondingly, 
assumed  uncorrelated,  zero  mean  white  Gaussian: 

E  [w  (k )  -wT(j)  ]  =  Q(k)6kj 

and 

E  [  v  (k )  *vT(j)  ]  =  R(k)6kj 

and 

E[w(k)]  =  0,  E  [  v  (k )  ]  =  0 

where  Q(k)  and  R(k)  are  covariances  of  excitation  and 
measurement  noise.  Now  if  x(k)  is  the  estimated  state 
value  after  the  kth  measurement  and  x(k|k-l)  is  the 
predicted  value  of  the  state  before  the  kth  measurement 
we  have : 
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x(k|k-l)  =  f (x (k-1 1 k-1 ) , k-1 )  ,  where  f  is  any  function. 
The  state  error  vector  is  defined  to  be 

x  (k )  =  x  (k )  -  x  (k) 

and  the  predicted  state  error  vector  is  defined  to  be 
x*(k|k-l)  =  x  (k  |  k-1 )  -  x  (k)  . 

The  covariance  of  state  error  matrix  is  defined  to  be 
P(k|k)  =  E [x (k) *xT  (k) ]  , 

and  the  predicted  covariance  of  state  error  is  defined  as 
P ( k | k-1)  =  E [x (k | k-1 ) • xT (k | k-1 ) ] . 

The  state  excitation  matrix  is  defined  by 
Q(k)  =  r(k)E[w(k) -wT(k)  ]  -rT(k)  , 
and  the  measurement  noise  covariance  matrix  is  defined  by 
R(k)  =  E  [  v  (k )  ■  v T  ( k )  ]  . 

The  Kalman  filter  equations  are: 

Pi  •*'  '  =  0  (k )  P  (k  |  k )  oT  (k )  +  Q  (k ) 

G(k)  =  P,  k-l)HT(k)  [H(k)P(kjk-l)HT(k)+R(k)  ]_1 
P  ( k  |  k  )  =  U  G(k)H(k)  ]  P  (k  |  k-1 ) 
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magnitude  of  the  gate  (-y/Z^  k-1)  +  R)  resulted  in  changing 
the  values  of  the  gains.  The  above  correction  makes  the 
filter  more  adaptive  now  and  the  errors  EB1  and  EB2  appear 
smaller  on  an  average.  The  representation  of  the  above 
correction  in  terms  of  trajectories,  gains,  predictions  and 
errors  appear  in  Figures  6,  7,  8,  9,  and  10  and  Table  4. 

The  computer  program  appears  in  Appendix  A. 


X  (1/0 )  = 


In  Figure  1  we  have  a  representation  of  the  best  track 
observed  and  the  K.F  track.  An  error  ellipsoid  at  the  25th 
stage  of  the  process  appears.  It  seems  that  K.F  follows 
the  observed  track  more  closely  than  the  best  track.  This 
is  due  to  the  adaptive  gating  and  Q  relative  to  R. 

The  gain  history  for  G(l,l)  fluctuates  between  0.7-0. 8, 
never  arriving  at  a  stable  value.  The  G(3,l)  reaches  a 
stable  value  of  0.12. 

The  above  appears  in  Figure  2.  The  track  for  prediction 
appears  in  Figure  3. 

The  errors  EB1  and  EOB1  represent  (YH-BLAT)  and 
(LAT-BLAT)  respectively.  They,  along  with  EB2  and  EOB2 
appear  in  Table  3  and  Figures  4  and  5. 

Implementing  Julian  Time,  we  have  a  comparison  in  common 
time  for  24  points  only.  It  can  be  seen  that  in  terms 
of  the  latitude  error  the  K.F  is  close  to  the  best  track 
values . 

In  an  attempt  for  better  performance  the  Q 
matrix  was  changed  in  the  algorithm.  Also  the  innovation 


errors,  in  terms  of  latitude  and  longitude,  exceeding  the 
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however,  does  not  have  its  major  and  minor  axis  aligned  with 
the  coordinate  system.  Instead  its  axis  (x1,  y'),  comes 
from  the  following  transformation: 

x'  =  xcos0  +  ysin0 

y'  =  ycostf  -  xsin$ 

where  x,  y  are  the  latitude  and  longitude  in  our  case  with 

9  =  1/2  tan-1 [  2|ov(x'y>] 

<7  -  a 

x  y 

where  cos(x,y)  =  p  (k  k-1) 

cr^  =  P11(k|k-1)  and  =  P22(k|k-1) 

The  new  variances  are  calculated  by 
2  2 

_»  2  _  _ Z  +  cov  (xy ) 

°x  2  sin2  6 

2  2 

,2  _  x  y  cov(xy) 

a  y  2  “  sin2  5 


Incorporating  the  above  equations,  error  ellipses  are 
presented  in  subsequent  figures  with  the  satellite  tracks. 


III.  ERROR  ELLIPSOIDS 

The  error  covariance  matrix  in  each  stage  of  a  Kalman 
filter  process  gives  insight  into  the  quality  of  the  track 
occuring . 

The  diagonal  terms  (P-q  and  P22^  are  t*ie  variances 
uncertainty  in  our  knowledge  of  latitude  and  longitude.  Their 
respective  off  diagonal  terms  are  covariance  between 
latitude  and  longitude. 

The  square  roots  of  the  diagonal  terms  gives  us  the 
rms  errors  in  out  estimates  of  longitude  and  latitude.  Having 
the  definition  of  the  structure  we  are  dealing  with  (in  our 
case  the  satellite  observations)  and  its  uncertainties 
(expressed  by  the  PCN  number-actual ly  by  the  values  of  R) 
we  can  see  how  the  K.F  performs  through  its  error  covariance 
matrix.  Expressing  the  P  matrix  in  an  ellipsoid  of  constant 
probability,  one  obtains  a  visual  appreciation  for  the 
worthiness  of  the  algorithm  parameters.  The  representation 
requires  that  the  errors  are  normally  distributed. 

The  joint  probability  density  function  is: 

e-l/2eT (kj  k-1) P_1  (k  j  k-1) e  (kj  k-1) 

where  e(kjk-l)  is  the  predicted  state  error  vector.  Setting 
the  exponent  equal  to  a  constant  value,  we  are  going 
to  have  a  curve  which  is  an  ellipse.  This  ellipse, 
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The  initial  condition  of  P  (error  covariance  matrix)  and  the 
Q  and  R  matrices  are  the  determining  factors  in  the  filter 
structure.  For  Q  having  main  diagonal  values  greater  than 
R  means  that  we  have  greater  uncertainty  in  the  state 
estimate  than  in  the  observation.  Thus  the  new  state 
estimate  is  more  dependent  upon  the  new  measurement  and  less 
related  to  prior  estimates.  The  inverse  is  also  true.  For 
R  having  greater  diagonal  terms  indicates  that  the  new 
measurement  are  subjected  to  stronger  corruptive  noises, 
and  so  should  be  weighted  less  by  the  filter.  The  gains  (G) 
are  lower.  The  P  is  responsible  for  the  initial  transient 
performance  of  the  filter. 


igure  8  K.F  Track  and  Prediction  Ahead  in  Latitude 
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Figure  9  Latitude  Errors 


Figure  10  Longitude  Errors 


V.  RANDOM  TRACKING 


To  examine  the  adaptability  of  the  K.F,  another  storm  was 
created.  The  equations  for  this  simulated  storm  were: 

BLAT ( s )  =  BLAT ( s - 1 ) -V  * T 
BLONG(s)  =  BLONG (s-1) -V  -T 

LAT(S)  =  BLAT(s)  +  V(s) 

LONG ( s )  =  BLONG ( s )  +  V(s) 

where  T=6hr,  V  =10°/24hr,  V  =5°/24hr  and  V  =  measurement 
x  y 

noise  (created  by  a  random  generator  subroutine) . 

Implementing  the  above  equations,  a  "true"  and  an 
"observed"  track  were  created.  These  two  data  files  were 
named  BESTRACK  and  OBSERVED  and  appear  in  Table  5  and  6. 

Simulating  with  the  above  new  data,  the  K.F  algorithm 
appeared  to  track  well.  Figure  11  shows  the  "true", 
"observed"  and  K.F  tracks. 

Figure  12  indicates  the  gain  history  in  terms  of  Gil 
and  G13.  This  approach  gave  stable  values  of  0.27  and 
0.01  respectively  after  the  4th  discrete  point  in  time. 

After  this  time  the  gain  does  not  vary  any  more.  This 
means  that  the  innovation  error,  (z  (k) -x  (kS k-1) )  is 


weighted  each  time  by  the  same  quantity  after  the  4th 
observation.  Having  lower  values  in  the  diagronal  terms  of 
the  Q  matrix  in  comparison  with  R,  in  this  case,  means  that 
we  have  greater  uncertainty  in  the  measurements  (observed 
data)  relative  to  the  model  uncertainties.  So  the  gains  are 
smaller  and  the  filter  no  longer  "tracks"  the  measurements 
closely . 

As  far  as  the  latitude  and  longitude  errors  are  concerned 
it  can  be  seen  that  EB1  (YH-BLAT)  and  EB2 (XH-BLONG)  are 
nearly  zero  and  are  smaller  in  comparison  with  E0B1  (LAT-BLAT) , 
E0B2  (LONG-BLONG) .  This  shows  the  ability  of  the  algorithm 
to  follow  the  "true"  values  more  than  the  measurement  if  the 
latter  has  been  corrupted  with  noise. 

The  error  values  and  the  plots  appear  in  Table  7  and 
Figures  13  and  14. 
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TABLE  5 


► 


BESTRACK  DATA 


TIME 

8L  AT 

BLONG 

1300. 00 

160.00 

3.80 

1806.00 

157.50 

2.30 

1812.00 

155.00 

1.80 

1813.00 

152.50 

0.80 

1900.00 

150.00 

-0.20 

1906. 00 

1 *7.50 

-1 . 20 

1912.00 

1*5.00 

-2.20 

1913. 00 

1 *2.50 

-3.  20 

2000. 00 

1*0. 00 

-4.20 

2006.00 

137.50 

-5.20 

2012.00 

1 35.00 

-6.20 

2013.00 

l 32.50 

-7.20 

2100.00 

130.00 

-8.20 

2106.00 

127.50 

-9.2  0 

2112.00 

125.00 

-10.20 

2113.00 

122.50 

-11.20 

2200.00 

120.00 

-12.20 

2206.00 

117.50 

-13.20 

2212.00 

115.00 

—  1  *  •  2  0 

2213.00 

1  1  2  .50 

-15.20 

2300.00 

110.00 

-16.20 

2305.00 

107. SO 

-  17.20 

2312.00 

105.00 

-18.20 

2318.00 

102.5 a 

-19.20 

2*00.00 

100.00 

-20.20 

2*06. 00 

97.50 

-21 .20 

2*12.00 

95.00 

-22.20 

2*13.00 

92.50 

-23.20 

2500.00 

90.00 

-2*. 20 

2506.00 

37.50 

-25.20 

2312.00 

85.00 

-26.20 

2518.00 

82.50 

-27.20 

2600.00 

80  .00 

-29.20 

2606.00 

77.50 

-29.20 

2512.00 

75 .00 

-30.20 

2613.00 

72.50 

-31 .20 

2700.00 

70.00 

-32.20 

2706.00 

67.50 

-33.20 

2712.00 

65.00 

-3*.  2  0 

2713.00 

62.50 

-35 .20 

2eoo. 00 

60 .00 

-36  .20 

2906.00 

57.50 

-37.20 

2312.00 

55.00 

-38 . 2Q 

2818.00 

52.50 

-39.20 

2900.00 

50 . 00 

-*0.20 

2906.00 

*7.50 

-* l .20 

2912.00 

*5.00 

-*2.20 

291 3.00 

*2.50 

-*3.20 

3000.00 

*0 . 00 

-**.20 

3006.00 

37.50 

-*5.20 
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TABLE  6 


OBSERVED  DATA 


ti*e 

1800. 00 
IB06. 00 
18 12.00 
1819.00 
1900.00 
1906.00 
19 12 .00 
1918.00 
2000.00 
2006.00 
20 12 .00 
2018.00 
21 00.00 
2106.00 
21 12.00 
21 13.00 
2200.00 
2206.00 
2212.00 
2213.00 
2300.00 
2306.00 
2312.00 
2313. 00 
2*00.00 
2*06. 00 
2*12.00 
2*18.00 
2500.00 
2506-00 
25  1 «:  .  00 
2518.00 
2600. 00 
2606.00 
2612.00 
2613.00 
2700.00 
2706.00 
2712.00 
2713.00 
2800.00 
2305 .00 
2912.00 
2319.00 
2900.00 
2905.00 
2912.00 
2913.00 
3000. 00 
3006 • 00 


LAT 


160.00 

157.30 
15*.  92 
152.71 
150.03 
1*7.63 
1*4.95 
1*2.61 
139 .99 
137.26 
134.62 

1 32.39 

130.31 

127.37 
125 .1  9 
122.4* 
120. 19 

11 7.40 
11*  .3  7 
1 12.92 
110.03 
107.65 
10*  .33 
102 .96 

99.90 
97.65 
95  .29 
92.** 
90.5  1 
37.1  a 
a*. a* 

32  .40 
30.00 

77.78 
75  .08 
72  .50 
69.96 
67.22 

64.79 
62.75 
59  .8  3 

57.37 
5*  .  72 
52.43 
50.06 

47.91 

44. 38 
42  .58 
40.14 
37  .2  6 


LONG 

3.80 
2.60 
1 .72 
1 .0  l 
-0.17 
-1  .07 
-2.25 
-3.09 
-4.21 
—  5 . 44 
-6.53 
-7.31 
-7.39 
-9.33 
-10.01 
-11.26 
-12.01 
-13.30 
-14.33 
-14.78 
-16.17 
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-13.37 
-13.74 
-20.30 
-21 .05 
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-25.52 
-26.36 
-27.30 
-20.20 
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-30. 12 
-31  .20 
-32.24 
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—34. 4 1 
-34.95 
-36.37 
-37.  1  3 
-33. 48 
-39.27 
-40. I* 
-40.79 
-42. 32 
-43.1 2 

-44 . 06 
—45 . 44 
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Figure  11  Trajectory  of  Fictitious  Storm 
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TABLE  7 


ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 
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7  j  tl  7  -  7 

0.17 

0.1  7 

:  .  4  1, 

0 . 4  0 

7  _  0  /-.  r. 

0.11 

0  .  1  1 

-0  .  1  0 

-  0  .  10 

7 ,  -)  7  7-, 

0  .  1  2 

0.A 

0  .  1 

0.  1  > 

7  2  >  7  i) 

0 .  1 

0 .  1  •; 

.  .  ••') 

) .  7o 

7  >  7  >»t, 

0.12 

0.12 

-  1  .  0 

-  0 . 0,J 

7  j  ■)  70  2 

0.27 

'» .  ’  J 

)  .  A  1 

0.01 

A  l  74,; 

0  .  10 

n  .  1  0 

-  j  .  1  > 

-  7  .  7? 

'A'i  on  4 

0 . 0  ? 

0.02 

-0.1 

-0.1  7 

7  _•  n  j  i  u 

-0.02 

-0.02 

-  J  .  1  0 

-  0  .  1  ,) 

7-OMln 

-0.07 

-  0 . 0  J 

0.0 

0 . 0 

AlfV,' 

0 . 0  4 

0 . 0  4 

' 

0  .  ->  H 

7,.  >  02  ■! 

0 .  or 

0  .  )  ■ , 

0  .  0  * 

0.0. 

7  2'  )  <  4 

0.04 

n.N 

;■  .0 

0 . 0 

7_  0  1.0 

0  .  O  1 

0.01 

-  0  .04 

-  0  .  O  4 

AM'm 

-0.07 

-  0  .  V7 

-  '  .  2  7 

-0  .  ;>  i 

7  J  >0,2 

-  0  .  1  ? 

-0.12 

~7.il 

-  .) . ,?  1 

7  „>■>  ,,  ’■  4 

-0.04 

-  0  .  r,  4 

.  .  2A 

-7  .  7  » 

7  —  1  rt  *'  4 

-  0  .  0  u 

-  O  .  .)o 

-J  .  1  7 

-0.17 

A  0„.  70 

-  (j  .  0  0 

-0.0  , 

0  •  0  7 

0 . 0  7 

7_>  >  >7r, 

-  0  .  1  ? 

-0.12 

-  )  .  2-’ 

.  ;>  1 

/*_  J  O  -  .J 

-0.  1  l 

-  0  .  1  l 

-0  .  <>  7 

-0.07 

7_  0  4.  *  >1 

-0.00 

-O  .  (1,; 

0  •  0 

0.O, 

A'<  1  >■. 

0 . 0  r" 

0.0), 

I.'  .  4  1 

0  .  4  1 

7  c'-  >7,; 

0 . 0  2 

0.02 

-4  .  1  4 

-O.I2 

7  .'■) 

J  .  04 

0.0  4 

0  .  0  M 

O  .  0  1 

7^0712 

0  .  0  1 

O  .  7  1 

w  •  1  4 

0  .  1  4 

A'lWI  4 

-  0.  1  -< 

-  0  .  I  -1 

—  •  *-  4 

-0.2. 
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J=  1 ,1) 


00  --  )  1  -  1  ,  L 


♦  i x,  ».<al  pfcgpav) 

lui)  ,1=1,5) 

1  -*  il  A  T  (  r>  .(  ,  '::ULr.1  I  3  £  N  I  I  ■;  1  C  A 


-  t  r  .  •  »  ^  v 


A  a  J.  '  .  . 


,57,  -jU) 


A  ’  1 1  £  (  J  ,  1  4  J  ) 

1+3  rCr 1AT ( *  0 • ,7  3  (  • * • )  ) 

•v"-  :£  (5,33) 

DC  1111  I  =  1 ,  N 

iiii  «Pir£  {b,  31)  (?n:  (1,3)  ,j=i,n) 

-5ITE  (3,  -t  1C) 

W  F  1 1 E  (  3  ,  1  5  0 ) 

1  50  FUP.-1  AT  (5  ( ,  '  ENT  ZF  1  HE  DI.1ZNSI0H  Cc  F  A  N  D  C  3  INPUT 
1  VZCICr  (i)') 

FEAD  (5,+Q)  L 
«  H  I T  E  ( b  ,  1  ■)  C ) 

loO  FU-MAT  (as,  »£‘HEa  I  hi  ELEVEN  IS  JF  TJE  PISTE  IE  J 1 1 C  X  '  , 

+  ix,  •  :i  \  1 3 1  :>  — I  a  li  ’i  a  . ') 
c :  -i  i  =  i ,  n 

D  C  4  J  =  1  ,  L 
a  F  I  I  E  ( S  ,  17  0)  1  ,  J 

1  7  G  r  J  ?  -•!  a  i  ( 5  <  ,  '  3  A  *  .V,  ( •  ,  I  1  ,  •  ,  •  ,  I  1  ,  • )  =  • ) 

P  E  A  D  (  5  ,  7  J )  I  E  L  ( I ,  J ) 

4  CONTINUE 
*31  IE  (d,  13C) 

UG  FCP./.A!  I*  J'  ,‘ja,  'THE  J  A  ii  *  A  SaTPIX  ( 3  I JT  P  I  3  1 1 1 C  N 
1  «Ai?IX)  •) 

EG  b  1=1,  M 

5  .  I T  F.  (  t  ,  1  J  )  (PEL(I,J),J=1,L) 

me  ■* p i  r ~  ( d ,  i ) o ) 

3  10  A  E  A  3  ( 5  ,  1  1  3 )  I  A  A 

i.7  (lA’l.^p.IT)  >  T  r  2 
I F  ( -  a  N  •  :i  z  m  j.  Y )  u  G  1 C  11  J  3 

*  '  1  1  -  (  "5  ,  1  “  "  ) 

r  E  A  D  (5,  1JI)  I  ,  J 

*  i  L  »  (  V  ,  1  7  J  )  A  , 

*’  r,  A.  ^  (  J  ,  *  'j )  *-  E  I.  (  l  ,  J  ) 


49 


DC  2  I =  1  ,  N 

2  'vEviTE  (a,?Oj  (  F  h  I  ( I ,  J  )  ,J=1,X) 

90  F  G  it  J  A  T  <  1  9  7  £  1  1.4) 

91  F  C  9  .  -I  A  T  (  1  P  7  Z  1  1 .  4 ) 


C  GO  IQ  683 

10  00  WRITE  (J,  1 0 G ) 

100  FCSXAT  (5  C,  '  CO  YCJ  WANT  TO  CHANGE  ANY  2 

1  »  A I  a  I  X  ?  ' ) 

900  ?E AD  (5,  1  10) IAN 

110  FC  F  M  AT  (A3) 

IF  (TAN.Zg. IZ) GCTu  19 
IF  (IAN. NS. IY) GOTO  1JOO 
n  ?  1 1 2  { d  ,  1  c  C ) 

1  z  9  I  OF '1  AT  (b  X  /  ' HI  CH  ELEMENT  OF  THE  MATRIX 
1  TO  CHANGE? 1 ,  / , 

OX,  *£J!1E1  AS  IJ;  KHEE5  I  IS  IHS  ivCw  AN 
1  COL'J'IN.  •) 

PE  AD  (5,  130)  I,  J 

130  r  Q  a  0  A I  (211) 

■< Pill  {3,-0)  1,0 
F  E  A  0  (  3 , 7  0 )  ~  E  I  ( I ,  J ) 

WRITE  (6,30) 

CO  i  1=1,  N 

3  h  PITS  {3,  30)  (PH  I  (I ,  J)  ,J=  1 ,  J.  ) 

WHITE  ( o ,  140) 

140  EG  Ml  AT  (5 X ,  'ANY  CTiiZP  CHANGES?') 

GCTC  9  JO 

14  ivPIIE  (9,  141) 

141  rOZ.IAT  (’  1  '  ,3X,  '  DISCRETE  TINE', 


C  F HAD (5,4423)  I 

C4  4  2  3  F CPU  AT  (1 1) 

L  =  1 

Vi  RITE  ('2,7731) 

7751  F  G  ?.  !1  A I  ( 5  X  ,  '  E.MSF  .K£  5  .MATRIX  • ) 

DC  7752  I =  1 ,  N 

DC  7732  J=  1 ,  L 
WRITE  (5,775  3) I ,  J 

7  75  3  FCF3  \T  (5X, » E  (• ,1 1, • , • ,11, • ) =  ') 

READ  (5, 7  77  c)  c  (I,  J) 

7752  CCNIIMtJE 

•«  R 1  1  7  (9,  7754) 

Vi  51  I*  (9,  77  :4) 

7  73'+  FUi'.l  AI  (57,  '  ICE  E  MAI 31  X  ’) 

DC  773  5  1=  1  ,  N 

WRITE  (9,97)  { j  (1,3)  ,3=1  ,L) 

7  78  5  v.  5 1 V  E  ( J ,  )  0 )  (E(I,J),.J=1,L) 

Co  <n id 

15  «  R 1 T  l  ( 3 ,  5  u ) 

oO  FORMAT  (5X,  'ENTER  I  HE  cLE.1E.iT3  DE  THE  IRAfiCII 
1  71  A  E  i<  1  X  -  -  P  H I  ’ ) 

DC  1  I =  1  ,  N 
DC  1  J  =  1  ,  N 
it  j  i  i  i  ( d  ,  'j  ) )  1,3 

CO  FO'H.AT  (5X,  »  till  ('  ,  I  1  ,  ’  ,  '  ,11  ,  ’)  =  •  ) 

F  E  A  l>  ( 5 , 7  0 )  F  h  1  ( I  ,  J  ) 

7  u  F  C  r 


N  D  =  1 


DC  1838  1=7,.'! 

Z  (I)  =0. 

EC  7  d  3  9  J=1,N 
hi  (I, J) =0  . 

HI  (I ,  I)  -  1  • 

?KI  (I, J) =0. 

Pal  (1,1)  =1  . 

K (I, J) =0. 

E (I, J) =0 
7b  j3  A  (I,  J)  =  0 . 

«  r  i  r  ^  ( .'5 , 4 10) 

WRITE  (8,777  1) 

77  7  1  FORMAT  (5  X ,  •  EC  YO’J  V.  A  N  T  TC  COMPUTE 
1  LINE  F?C2!  '  , 

+  1 1 ,  '  A  S  3  ?  • ) 
h  E  A  D  ( 5 , 7  7  7  < )  T  A  N 
117  2  FC?.1AT(U) 

IF  (IAN  . E*. IPH) GC  TC  7  7  73 
IF  (IAN  .  NE.I  Fii)  GO  TC  1  8 
7  7  7  J  WRIT 2(8, 7774) 

77  7h  FOhdAI  (//,5T,  ' EST2T  THE  A  MATRIX') 
DO  7775  1=1, N 
DC  7  775  J  =  1 , N 
*  F i I E  (8, lilt) 1 ,  J 

777b  FORMAT  (ja,  • A  ('  ,11, ' ,  1 ,11,')  =') 

F  E  AD  (5,7778)  A  il ,  J) 

111  6  t  G  !u1  A I  ( F  1  0 . 3  ) 

7775  CONTINUE 

WRITE  (9,  7771) 
hi?  I  IE  (d  ,  77  79) 

77/9  F  C r  1  A  T  (  5  X  ,  *  THE  A  MATRIX  ') 

DO  1 1  a  <j  1=1..! 
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•,  SITE  (8, 27) 

r  C  3 .1 A I  t/2X  ,  1  '/, 

*2X,  'DO  NCI  USE  E  A  1  KEY  1C  EXIT  ?F:GGSA1  UNLESS 
1  Y  C  J  W AMI'/, 

* 2 X ' T 0  LCSE  YCUL  INPUT  FILE  I HAT  IS  PRINTED  ON 
1  UNIT  <4.  '/ , 

#i *******#*******************************$£*******1 / f 

*2  X  *  FOR  CHECK  CUT  ACTIVATE  THE  GO  IC  d  cJ  3  STATE M551 

1  CSV. 

*  2  X  'LINE  135  TC  END  If  17  HC  01  iBc.OP'/# 

*  2  X ,  • DO  YOU  NANI  TO  DEVELOP  THE  INPUT  FILE  ( Y/N)  7 ' ) 


READ  (Jf  2  3) IANS 
33  FORMAT  (A4) 


IF  (IA..S  .  iNE.  IUSt:  1 )  30  10  3999 


v. 

c 

C 


77  77 

FCSM. a  I  ( 1 H  1 ) 

WRITE  ( 6 , 

1224) 

1  I  j  U 

FORMAT  (  ' 

ALL  INFUL 

WRITE  (  3  , 

U35) 

1  2  j5 

EOF..'!  A  I  (  ' 

1.  HERE  CT 

1016 

FORMAT  {1 

H  0 , 5  X  ,  5  H 

19  32 

WRITE  ( d , 

13) 

1  3 

FC'J  IhT  (/ 

,  5  X  ,  '  THE 

WELLE  ( 8  , 

20) 

JO 

FORMAT  (5X,  »  E 

1  (UP  TC 

3)  .*) 

ENT  EE  THE  C*DER  OF  ThE  SYSTEM 


LEAD ( 5 , 4  J  )  N 
'4  0  F  J  K  i  A  I  ( I  1 ) 


_  ^ 
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I 


t_ 

C  t»KK  (I,  J)  =  i  (K/K)  ,  (CG  V  EFPO?  AT  K  -ilV 2Ti  X  SA/PLES) 

C 

c 

C  PKKK1(I,.7)  =  ?  (K/K-1)  ,  (CCV  EFPO&  AT  K  GIVES  K- 1 

CSAJ5ELZS) 

C  N  =  N rJ  13  3?  OF  P 0 VI 5 ,  M  =  N'JSBE*  CF  CCL.,CGS. 

C,NC  A  :«D  .1C  A  E! 

C  !«G  =  NU16Z9  OF  TTECAUCNJ  Cr  THE  SILTED 

C 

C 

(J  *********  i<F'n'  NEnS  ************************************ 
WRITE  (6, 7  17  1) 

7171  FCE.'iAi  (/,2X, **************  HEW  NEWS**************'/, 

*  2  X  ,  '  T  H  Z  oTCF.,1  HACK  INPUT  VALUE’S  Alt  E  A  VI AL A 3L2  '  /  , 
*2X,'THES2  Air  E  SC:1C  PBINTE2  TC  THE  TERHIi.AI  AT  TEE'/, 
*JX# '233INN1KG  CF  THE  PFCGSA.M  —  ALSO  SEE  LISTING'/, 
*2X,  'IF  THE  t-itCGRAi';  ENDS  NO?  M  ALLY  AN  I«FUI  FILE 

1  HI  LI'/, 

*  2  X ,  ' 3  E  PRODUCED  T H A  *  HA 1  PE  PRINTED  CF  USED 
1  FOR  INPUT'/, 

*  c  X  ,  'USING  1 h E  SAMS  FORMAT  STATEMENTS  IC  READ  AS 
1  WERE'/, 

*  2  X ,  '  USED  TC  »F  ITS  CN  UM1  4.  SC  WITH  A  FEW 
1  HCCIr  ICAlICN'b  '/, 

*  IX,  'THE  cP ARCH  AFCUNO  IIiE  INPUT  CAN  HE  USED. 

1  PPESERILY  '/, 

*2X ,  '  A  PILE  CF  THE  INPUT  DATA  IS  USING  ’’FODUCZD  AND 
1  WILL'/, 

♦2X,'.jE  FOUND  AS  --  K  CUIPJI  A--  CN  YCUR  A  -  CxSX.') 

C  ****  TiiS  FIRST  CUESTICN  TC  THE  T  E  ~  *1 1 N  A I  — ASKS  IF 
C  AN  INPUT  FILE  IS  IC  lE  USED 
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J  7  3  ESI. AT  ( 2 X  ,  7  F  1  0 . 2 ) 

C  •■FITE  (  J  ,  -i  7  )  (TINS  (I)  ,  L  A  1  (I)  ,LONG  (I)  ,  ETUI  (I) 

1  ,  B  L  A  I  ( I )  , 

C  *  3LC:,u  (I)  ,  d  X  I N  D  (I)  ,  1=  1  ,  M  2) 

C  V  HITE  (  i  ,  3  7  4)  (I  IMS  (I)  ,LA1  (I)  ,LGa3  (I)  ,  EIIM?  ( 1) 

C  1  ,  2L  AT  (I)  , 

C  *  PC N  { I )  ,  1  =  1,  NZ) 

C  J  7  4  F  0  3  A  A I  ( 2  X  ,  6  E  1  0  .  2 ) 

C 

C  THIS  PROGRAM  COMPUTES  TEE  FOLLOWING  KALMAN  PILISr 

C 

C  -  1 

C  G  (K)  =  £  (K/K-1)  *hl*  (H*P  (K/K-  1)  *HT«-F) 

C 

s# 

C  ?  (K/K)  =  (I-G(K)  *H)  *F  (K/K-1) 

C 

C 

C  P  (K/K-1)  =  F Hi*?  (K/K-1)  *PH11+C 

C 

C 

C  g  (I, J)  DEFINES  THE  COVARIANCE  OF  THE  PER  SAMPLE  RANDOM 
C  EXCITATION  Or  IKE  P5CCE53 

C 
C 

C  Ef  (I,  J)  DEFINES  THE  PAN  DOM  (GAUSSIAN)  M  E  A  S  ‘J  t  Z  .1  ENT  NOISE 
C  WHICfa  IS  ADDED  IC  THE  C SS£ 5 V  A  cl E  SIGNALS 

C 

c 

C  { I  ,  J )  IS  THE  I  u  £  tN  i  I T  i’  MATRIX 

r 

C 

C  1 1  =K  THE  DISCRETE  POINT  IN  I  IMF,  THE  STAGE  Or  THE  Di'CCrSS 
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*'»  i.  ",  O  J  -> 


i  y  i  n’  i 


/  •  Y  •  / 


(j  ^  7  =  NC.  Of  C  2  5  £  E>  V  L  J  V  A I  J  -  3  (£.\Ij-L-.E) 
C  M  Z  =  A  C  .  Or  E  2  i  I  I  3  A  E  VALUES 


MZ  =  53 
JZ=50 

DC  132  1=1,4 
F  K  K  ( I ,  I )  =  1C  0  0. 
13  2  ?kk:ii  (I, I)  =1000. 

{  1 ,  1 )  =  .  0  0  0  0  0  1 
,»  (2,2)  =  .  CGCOO  1 
C  BEAD  OBSERVED  VALUES 


L 

c 


BEAD  (2,  1  1)  (TIME  (I)  ,  LAI  (I)  ,  LONG  (I)  ,PCN  (I) 
11  FOB MAT  (7a, 5 6. 2, F3.1, IX, 54.1, IX, £1.0} 

F  E  AD  (2 , 1  1 )  (TIM  2(1)  ,  L  A I  (I), LG. '13  (I)  ,FCN  (I) 


1=1, NZ) 
1=1, LZ) 


C 

c  LEAD  j£SI  I  FACE  VALUES 
C 

p  E  AD  (3,1*)  (ETIME  ( J)  ,  2  LAI  (J)  ,  BLONG  (J)  ,  c  *  1  N  D  (J) 

1  ,  J  =  1  ,  M  2 ) 

14  FCrMAl  (oX,F4.0,F4.1,1.(,F4.1,2X,F3.0) 

C  BEAD  (3,11)  (LIIME  (J)  ,ET  AT  (.1)  ,  BLOND  (J)  ,?CN  (J)  ,0=1 

C  11  FC.-.MAI  (<.X  ,4510.2) 

C  11  r  C AT  (-iX  ,  f  *-  0 , 54  .  1  ,  1  X  ,  t  4  .  1 , 2.< ,  E  3. 0) 


C  ECHO  VALUES  *  *  *  *  Z  EMC  7  E  HI  XT  TliBsE  LIN^S  TO  ELi/K-A 
C  EC.iT)  °  B I N  I  ChfcCK  *********** 


KBITS  (0,373)  (TI.1E  ( I )  ,LAI(I)  , L3  S  i  ( I )  , E  T 1 M  E  (  T ) 

1  f3LA«.(l)r 

♦  3 LONG  (I)  ,  rCN  (I)  #  1=1,  N  3) 

I.  B  I  IE  (d,  3  7  j)  (TIME  ( I )  ,  L  AT  (I)  ,  LONG  ( I)  ,  PCN  U  )  •  1=  1 
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APPENDIX  A 
COMPUTER  ALGORITHM 

KALIAN’  rILTZ." 

DI 1ENSICN  HI (12, 12)  , H (12,12), R(12,  12)  ,  G  (  12,  12) 

1  ,  ?iiIZ  (12,  12)  , 

*  0(12,12)  ,  G  3  1  (120)  ,criK(12,12)  ,  ?  K  X  >1 1  (12,  1  _  ) 

1  ,  3  11  (  U,3)  ,G22  (  123) 

DI  -1  ENSIGN'  IFEAD(1'J),I»biITE(10),Y(12G),Y:i(120) 

1  , XH (120)  , G*c  {  12  0) 

DIMENSION  CEL  (12,  12)  ,  A  (1c,  12)  ,3  (12,  12)  ,D1  (12,  U) 

1  ,  D2 ( 12 , 12) 

Cl  MENS  ION  2  (12)  ,2  (12)  ,GE  (12)  ,  Xh?  (120)  ,  I!iP  (  120) 

1  (120)  ,  EX  (120) 

DIMENSION  DELI  ( 12,  12)  , PHI  (  12, 1 2), II  (120)  ,  X X  *  +  *  ( 1  ^ ) 


D  I  MEN’S  I  ON  JAM  (5)  ,D  (U,  12)  ,  XP  (25)  ,22(2  3)  ,  X  ILLS.;  (12J) 
(12,12)  ,?  (12,  12)  ,AI  (12,12)  ,  2 R !< -1 1  (12)  ,  IK  (I2J) 

1  ,  ?KX2 2  (  12  3)  ,  ? 22  (  1  20)  , 

*FT  ( 1  2,  1c)  ,  ET  (1  2,  1  2)  ,XNK  (  12)  ,  XKKM  1  (  1  c)  ,  ZXY  ( 12  3) 

1  ,  EO  3  1  (123) 

J  i.MENS  ION  TI.'.E  ( 123)  ,  LAI  ( 120)  ,  LONG  (12  3)  ,  FUME  (120) 

1  , 5  LA  I  (IzO)  , 

*  C  L  O  A  G  (12  3)  ,  E  ;*  I  N  L  (123)  ,  XX  '  (123)  ,  /  X  /.  >>  (  12  1 )  ,  “>*..(  (120) 

1  ,  £032  (1-0) 

INTEGER  I  X  X  7  2  (  1  0  3 )  ,  I K  '1 4  <3  (10  J )  ,1^2- t  (130)  , 

*  Ir,  X  12  (130)  ,  IK.1  j  (  1  00)  ,  IisT-o  (100)  ,  IK?  1  2  {  1  '3) 
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VI.  CONCLUSIONS 


The  K.F  approach  to  estimate  the  storm's  location  appears 
to  be  very  accurate.  That  comes  from  the  comparison  with 
the  meteorologist's  analysis  results.  Concerning  also  the 
fact  that  the  latter  was  performed  after  the  storm's  occurence 
one  can  see  the  advantages  of  the  K.F  algorithm.  During 
operation  of  the  filter  the  actual  residual  sequence 
[z  (k) -z,  (k|  k-1)  ]  is  compared  to  the  gate  [  \/P  (k  k-1)  +R  (k)  ] 
which  actually  is  the  square  root  of  the  residual  covariance. 
Being  a  white  Gaussian  sequence,  the  residual  is  bounded  by 
this  gate. 

When  the  gate  is  exceeded  the  model  is  determined  invalid 
within  the  filter  and  a  modification  in  terms  of  Q  and  G 
takes  place  to  adapt  to  the  situation.  At  this  point,  and  if 
the  excess  occurs  in  only  one  component  of  the  vector  residual 
process,  one  can  further  deduce  that  the  measuring  device 
generating  the  particular  component  is  the  source  of 
difficulty  (a  sensor  failure). 

The  error  ellipsoids  of  the  process  also  give  insight 
into  the  filter  performance  in  a  more  general  case,  referring 
to  many  sensors  with  a  greater  variety  of  uncertainties,  the 
adaptive  K.F  algorithm  could  be  a  very  advantageous  approach. 
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Figure  14  Longitude  Errors 


Figure  13  Latitude  Errors 
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•'"Ii't.  t  j,  3  19  1)  GAT2,  E  ( 1) 


9 1  y  i  !•  c :  \  i  (  1  e  f  r  c  ~ 

AT  GATE.  GATZ= 

’  ,  F  1 0 .  ♦ ,  9  X  ,  1 

'"(1; 

1  ,  r  1  0 . 4  ,  •  X  X  ’  ) 

'J  1  9  2  F  C  r  .1  A  I  ( 'J  X  ,  1  2  r.  E  0  E 

31  gate,  gate- 

» ,5  10.  4, ?X,  ' 

’  Z  (4) 

1  ,710. 4,  'XXX') 

C 


L 


V- 

G1 1 (X) 

=  J 

o. 

1) 

4. 

031 (K) 

=  O 

(3, 

1) 

c 

DC  3022 

1= 

1,3 

c 

»  7 1 T  r-  ( a , 

00 

)  (E 

X  K  Ml 

U 

,J)  ,0=1,3) 

3  J  22 

TRITE  (4, 

0  0 

)  (2 

K  X  A  1 

(I 

,J) ,0=1,3) 

2204 

••17  IT 

E  (9, 

59) 

’  '  TO  '  T 

E  (  0, 

5  3) 

J  0 

7  j  c  M 

AT  (/ 

1  a  a  :  ?  l  a 

vj 

1 

) 

DU  3 

02  4 

1=1,3 

'..TIT 

^  ^  • 

90)  ( 3 0) 

,0 

=  1, 

-i) 

302  4 

WH II 

E  (4, 

90)  (G (1,0) 

,  J 

=  1/ 

?) 

.HIT 

E  (o. 

3 0 )  (ZXXS  1  (0) 

,0  = 

1,0) 

Vi  71  ' 

E  (^, 

Eb  1 1 ) 

K  4  I  * 

E  £4, 

•-‘•11) 

a  3  11 

Fv.  3  .0 

AT  (/ 

1 

’) 

.  'UT 

E  {  i. 

30)  (2(1;  ,  0 

=  1 

,  :•! ) 

* .  i.  i 

E  (4, 

-  4  )  £  -  (0)  ,0 

=  1 

,  .v) 

C  A  T 

0  0 » . 

o  (G,  2,  j, 

1, 

-  t 

N  u  ,  1 D  ,  L  D ) 

c 

A  i  i 

4(4, 

(  G  4  £  0  )  , 

0  = 

1,0 

) 

J.iLL 

A  w»  -> 

i  x  ?  o  i ,  ;  z , 

A  , 

i,< 

K  K  /AD,  -ID) 

L 

J  <  ’i  i_  L 

?  1C 

C  ( t?  H  7, "KF, 

3  , 

i 

,  X K  K 3  1  ,  4  D 

rf  "•  1  A 

E  (  >  * 

4  0  11) 

»•  r  '  i.  i. 

E  ( 4 , 

c  oil) 
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r:  11 


x v ’ 


i 


c 


, :  :tl  (  i.  r.  2)  (a-'*  (j)  ,  j=  1,:.) 

(x.xk  (j)  ,j=i,n) 
in  r<)  =g  (i,  i) 

GjI  (K)  =1(3,1) 


?1  1  (A)  =  FXK:11  (1,1) 
P  K  1  1  (  K )  =  r  X  K  (  1  ,  1 ) 


Y  H  (K)  =  XFK  (  I) 
AH  U)=XKK  (2) 


KTIIMNl  (TIME  (_<)  ) 

2FTIH2=FLC  AI  (Kll^r) 

(PXTI-12.  Ji.71  IE  (?)  )  GOTO  3d  13 
( « )  =JUL;i!' 


L 

EO  13  1 

(K)  = 

=  Y  r. 

(K) 

-LAT  (X) 

C 

ZOO  2 

(■<)  : 

-- h 

(K) 

-LONG  (X) 

3d  13 

CO  NT 

I N  ’J ' 

C 

9PIT 

E  (9, 

,36 

n) 

C 

W  FIT 

E  (d  , 

,  3  3 

12) 

8  3  12 

F  C  4  1 A  T 

(/’ 

X  K  K  31  •  ) 

i.iirE  (o,  )0)  (Xkk-m  (J)  ,j=i 

PITS  (4,GU)  (A.xv;11  (J)  ,  J  =  1 ,  ;i ) 
Y  HP  (  F)  =  X 7  K  M  1  (1) 

X  H P ( K ) =  XKK3i  1  (2) 


c a  = 1 3  0  0 ♦  11*10  0 


•j n  5  DC  ?  o  8  h  1=1,4 

9  e  8  6  IP  ( 1 1 H  c  ( K  )  -  L  J  .  T.  a .  ) )  GC  I C  0  0  17 
tc=jjtc 


IP  (Ba-iI.IE(K).Li.O)  G  C  1C  ItMio 
0997  IF (32-11  IF  (K)  . Gl. 0)  GO  1C  9o99 
Y  K  li t,  (:<«*”)=  X K K  (  1 ) 


X  X  :i 6  ( K  'X  K )  =  XK«  U ) 

GC  Tu  v  ti  3  d 

:i*  n  cc  =  2--  (15  n*I  I*  1  1  1) 

A  1 1  I  ( 4 , 5  a  E  L! )  C  C ,  L.  r 
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j  i  ; 

j  1 


1 


;  {  l  '  *.  »  i  /  4  -j  •  i. ) 


'!  -■  l„ 


C 

095/0  04.1*2  J  J  T 1 .1  £  (IZ1AE,  JULtiF.) 

w  ****** ***************  i.oGi'ic  :■  c n 1 1 n z  ********** 

C  ***  CALC'JL  41  FS  JULIlM  T  I '•)  5  FPC.9  YEA?  190.) 

l  ***  i y r  =  y  t a :■. ,  ::'.c=«c4t:i  (ialch)  ,  ica=oay, 

c  *  *  *  i ;  • c  o  o  r  o  a  y 

I  Y  5  =  1 S  5  5 
1.10  =  5 

nA=ni  i£/  ioo 

:;!5=r:i,i  £-ida*  ioo 

CALL  9  'M  C  H  N  ( I  i  7  ,  1 .1 C  ,  I  C  A  ,  I  ri  5  ,  J  l)  L  ti  t. ) 

o :  .j  :j 

E  N  0 

.5'JiJRC>JT:  NS  N  n .1 C Z  ,N  (  1  SAP  ,  OC  ,  DA,  H"  ,  JOLH*  ) 

C  ***  CALI.  SC  lJY  S  J  E  /  C  rJ  1 1  N  E  JUIISZ 

C  ***  CALC  111  ALSO  JO  LI  AO  LAY  A  A  0  JCLIA.N  00  OP 

I  N  7  H  G  E  L  IN  IL  (1i)  ,  1 1  .1 ,  C  A  ,  :•  ? 

0  A  :  A  I N  I  0/  0  ,  5  1 , 0  -J  ,  v  C  ,  1  2  0  ,  1  5  1 ,  1  5  1  ,  4  1  2  ,  4  4  5  ,  2.  7  J  r  5  5  4 
1  ,554/ 

I E  =  ( Y  E  A  r  -  1 4  J  J . )  +3^5.25-0.25 

I A  D  0  =  0 

IE  ( IOC  (Y£  iH,  4)  .  GT  .  C)  GOTO  t>0  3 

I  A  Z  D  =  1 

».  J  j  J  0  L  0  A  =  1 S  I D  (  A  t ; )  C  A  I A  D  2 

I  ;!  ;•  =  5  4.  +  ( i  0  4  .]  rj  I  c  A  -  1 )  4  h  2  4  J .  5 

I  Y  t  A  E  =  Y  E  A  l  .  -  1  ■  j  ) 
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'  L  \  -  i  .1  iM-0~3  I..  ?  C  S  IT  I J  ’i  OF  .<  a  L  /.  A  .  FH.Ii:- 


DIC.  :uJ.£Av 


t  z.  J)  .  \  .  A  1  • 


/  A  *.  j  «  A  n  A  r  j  .  I  N  0 

z  a  t  a  i  c  HE  I-  I  L  £ 


TO  2=  1=  1 , 50 
1 3  T I  .*  !£  =  I  NT  (II  ^  i  (I)  ) 

.ALL  .miSS  ( 1 3  1 1  .<!  £  ,  J  0  I  ri  I-  ) 

1 3 J :j i  ( i)  =  j;j:  ht 

DO  IT  J  =  1 , 5  0 

tot  =  i;,t  (E-:i.i  *  (J) ) 

CALL  JUT  IMI-:  (IH5I,  UJ) 

I  3 •  J f J L  (J)  =  ILJ 

I?  ll&.ini  (J)  .  I5JUI  (I)  )  OU TO  13 
EL  I  (I)  =  Y  u  (I)  -  clAT  (J) 

Ec2  (I)  =  X  Li  (  L)  -bLC  NO  (J) 

TO  d  1  ( I )  =  Y  H  ( I )  -  L  A  T  ( I ) 

ZOz>  2  i  1)  =  X  H  (  1 )  ~  L  C  N  0  ( I ) 

EO  3  1  ( 1 )  =  L  A I  1 1 )  -  I  L  a  I  (.)) 
i\  Od2  (i)  =  L  ^  Ntj  (I)  —^LONG  (J) 

JTI’IE  (I)  =  IEJ  'JL  (J) 

CONILN-JE 
CO  NT  I  NOT 
••  ?:iE  (1  -*,2  3  3) 

A  £  1 1  £  (o,  290) 


F<i  £  1  AT  (//  ,  1  J  v  ,  Ml  I  2E  ’  ,  M  ,  •  2  2  1  •  ,  r;  X  ,  ’  £  r  2  •  ,  q  Y  , 

1  •  TOd  1  •  *  £CE2*  ) 

DO  313  ”=1,53 

IT  (J Ti  l-  (?)  3CTC  2  1  3 

(o,  2 j 3)  mtmf  (N)  ,  ?.2  i  (H) ,  122  (N)  ,  in_-i  (N) 

1  ,  £  0  3  z  (  J )  ) 

>■  I  r  £  (  1  *,313;  (Jr  ME  ( )  ,  E  3  1  (N)  ,M^(N)  ,  L  M  1  ('-) 


ACE  -.c.L!\  D  ( A , N, 1  ,  N  !),  ID,  I  A  i  u  ] 

31  v  PC.  3 10:.  A  (?iD,v!C)  ,11 E  AC  (10) 

Du  1)  I  =  1,N 

A  E  A  3  ( 5  ,  2.0)  (A  (1,  J)  ,  J  =  1,rt) 

F GP.1AI  ( 3  P  1  C  .  5) 

.1 E I J  E  :i 


A  OBkCUTiAE  1  Vi  n I  T  E  ( 4 ,  3  ,  1 ,  * J, 10, I  * « I T £) 

DUE  N3ICN  A  (  N  J  ,  1  C )  ,!*?::£  (10) 

30  io  i  =  i , 

v  y  r  e  [1,12)  ( : , ,  4  { i ,  j  ) ,  j  -  i ,  i ) 

'•*’  P  1 1  E  (3,  EO)  (1,0,  A  (I,J)  ,  j  =  I,') 

r.A/!  AT  [2  (JX,  •  (’  ,  II,  •  ,  •  ,13,  •)  =  •  ,  1  ?  E  1  0  .  3 1  ) 

.i  EE’n:, 


J  J  O  '  i  i  I  . 1  .*  .  J  "  (  J  L  H  *  ,  lO.  C  ,  i  i  /  |  C  -  1  -  £  ,  Aii  ,  (  ll  ,  i.  . 
i  LO  N  o ) 

L  !  E  V  3  :  •'  1  1 1  X  E  (1)0)  ,  J 1 1  X  E  (JO  3 )  ,  3  L  A  I  (JO))  ,  c  1. 3  K 

,  i :•  (  jco) 

[  -E;.3  io..  <H  (3  )  ))  ,  :  _•  1  (  j0  3)  ,  EE2  (303)  ,  IEJ'JL  (3  J  J) 
r  lA.I'II  (  VK)  , 

: : :  12  (  j  )■.' )  ,  so  l-  i  (  j^i  ,  zc.32  ( _> 0 d ) 

-  "  A  I.  *  4  L  AT  (j.OO)  ,  LENA  (  33  0) 


DO  li'j 

r 

X 

— 

1, 

‘i 

CO  1 20 

J 

= 

1. 

1 

o  { I  » J) 

= 

A 

u. 

J) 

h  E  I‘J  P  N 

IF  t;+1 

. 

1  5 

,14, 

00  1  4j 

I 

= 

1, 

a 

DO  1  40 

J 

= 

1  , 

'i 

0 (I, J) 

= 

- 

A  (I 

,J) 

Rirj?N 

DO  150 

I 

= 

1, 

a 

00  150 

J 

= 

1  , 

A 

C (I, J) 

= 

K 

*h  ( 

I,  J) 

?  S I  u  l-i  :i 

E  S  0 

SO  3  F  COT  I') 

E  ^  17  C I  P  (  N  ,  t  ?  ,  A  , 

DIMENSION 

A  {.-.,!!)  ,  X  (.1,0) 

DO  1  I  * 

I,:i 

uO  1  J  = 

i,:i 

x  £1,0)  = 

0  . 

0  0  o  »  = 

1  , 

= 

1 . 

DO  ja  L 

=  1,-1 

K?  =  0 

2  =  0. 

DO  U  K  = 

2,  H 

i (2.0c. 

AKS  (A  (K,  J)  )  )  oc 

l,  -  Ac3  (  A  (  K  ,  L)  ) 
.<P  =  K 


0  0  J  I  i  N  J  L 


C  2-4  a?1T:  (2,  iO)  ITE.'IFI  (I  ,0)  ,0=1  ,  1) 

c  a  _l  ...  d  u  (T  e/.  p  i ,  :• ,  .1 ,  .i ,  r  5  2  r  1 , d  ,  :i  d  ) 

CALL  ?5Ci?  tv»0.  OK  COO  1 , T 7 2 ?  1 , TEH  ?  2  ,  K  i :J. 
c  £21?2  (1,  (c,z)/DEC 

C  IE2<?2  U,  1)  =-IE2P3  U,  1) /DEI 

C  I7«?2  (1,2)  =-IE:i2J  (2,1) /CEL 

C  TE.1P2  U,  1)  =-TEri?3  ( 1,  1)/D£I 

0  WRITE  (3,31) 

c  DC  27  1=1,2 

C  z.  7  WRITE  (8,  30)  (TE3P2  (I  ,  J)  ,  2=1  ,  '.•) 

3  1  rCo  1  AT  (•  (  H  ?  H  +PJ  -  1  ' ) 

17  (KEE-2)  101,113,  10  1 

110  A' PILE  (3,111) 

111  FORMAT  ( 5  ii  K  E  F  =  2 ) 

101  CALL  ?  7  C  D  ( '1 E  2  P  ,  T  E ! !  P  2 , 8 , 2 , 2  ,  G  ,  3  J  ,  M  D  ,  L  D ) 
C  ii C  I  E  HEFS  P i\ K  ( I ,  J )  =  7  (K/K)  3  HERE  ?(K/K) 

C  (I-J  (X)  *H)  *P  { «/K-  1) 

wALL  ??CD  (G,h,  :;,.1,'I,:£2P,ND,  «E,IC) 

C  WRITE  (3  ,  30  ) 

C  30  FORMAT  {'  G Ii  '  ) 

C  DC  25  1=1,  N 

C  2  8  «  3 IT  E  ( 3 , 0 0  )  (I2M?(L,J)  ,  J  =  1 ,  N ) 

DC  10c]  I  =  1  ,  >1 
DC  1)8  J  =  1 ,  N 
1  Jo  1  E  2  (I ,  J)  =  -T  F.  1?  (I,  J) 

C  WFITE  (o,  37) 

j7  FORMAT  ( '  HI  ») 

C  CG  45  1=1,3 

0  4  3  WRITE  [6,  K)  (  21  {I,  J)  ,0=1,3) 

I,Tc N, N , TIME , N J ,  ID) 


C  A L  L  ADD  ( fl 


1  ,  .'I  J  ,  L  D  ,  v  ! 

c :  -i  z'iz  ir' »  u,  U)  (1^,  12) (12, 12) ,  j  (12,12) 

1  ,-(12,1^), 

1  .ilv12,12)  ,H2  (12,12)  ,  IF  *?(  12,1^)  ,IE.1-2  1(12,  U)  , 

i  (12, 12) , 

KPHi  (12,12)  ,  FHIT  (12,  12)  ,  P  '<!<  11  (12,12) 

3  ( rT )  =  ?  ( K  /  K  -  1  )  *di*  (2*2  ( X/  K-  1)  *  H  I  +  S) 


?:(!*?  (A-  1/A-  1)  *LHII  +  C 
IF  (K.  S'2.  1)  CC  rc  sJoi 

CALL  TUANS  (?  !i  I ,  N  ,  N  ,  P  H  1 1 ,  MU,  MD) 

CALL  P.tCC  (?KK,  '‘HIT,  N,N,  N,  IFHP,  NC,V,C,  IP) 
CALL  ?3CC  (r  HI,  IF  HP,  N,M,  M,  12.1?  1  ,  Ml  ,HD,  LD) 
CALL  ADD  (IFXPl,  2.  ,  h  ,  N  ,  2  KK  A  1  ,  MD  ,  .ID) 

A2ITZ  (2  ,  t:  5  5 ) 


pd5  FORM  A'i  (/  '  .1 A 1 2  1 X  ZKXMl  ') 

DC  j)  i2  1=1,  M 

v?:it  (o,  2  ))  (?k:cii  (i,j)  , j  =  i , M) 
j  0  j  2  *FI2E  (4,  JO  j  ( FK  K X  1  (I ,  J)  ,  J=  1  ,  N) 
j  A  *j  L  C  C  2  1 .1  J  jj 

CALL  T  A  A  N J  ( h , M , N , HI ,  A D  ,  X  D ) 
c  i.  :-i:z  (?,  j*j) 

J3  FCCXAie  il  ') 

C  DC  2  2  1=1,  « 

C  22  WRITE  (6,'iO)  (H  (I,J)  ,-J=1,S) 
jO  FC?H  AT  ( 1  P  7  E  1  1.  4) 

C  „  RITE  (3,3c) 

jo  Zo  ESA  I  (  1  Li  I  '  ) 

C  DC  z3  1=1,  S 

C  23  «*  f  1 2 Z  ( b  ,  ) 0 )  (hi  ( i  ,  J )  ,  J  =  1  ,  X ) 

CALL  ?  jC  D  ( ?  <  K  :•!  1  ,  H I  ,  M  ,  N  ,  .M  ,  I  E  X  ?  ,  N  Z  ,  r 
CALL  2  ?.  0  0  (  H  ,  1  £  X  2  ,  M  ,  N  ,  1 , 1  E  A  P  1  ,  N 1 ,  X  D 
C  WHITE  (3, 3o) 

3  3  F  C  L  *  A I  (  '  1:  7  HI  '  ) 

C  CC  2  4  I  =  1  ,X 


,LC) 

LD) 
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=  :/Z. co*r;>: 


U  1  1 

Li  - 

1,3 

0  1  1 

I C  = 

1,-<- 

DP  (IF 

,IC) 

=  T/F*C(I2,n) 

HI  (IF 

,  IC) 

=  ?HI(IE,IC)+CC?.(I5,IC) 

fill  (I 

?,IC) 

=  rill  ( I E  ,  I C )  *■!/  (  (F+  1. 

*CCF  (1 F,TC) 

1  fir  1  (IF,  1C)  =  TEFC  (IP,IC)  +  1/  (^+1.  }*CCF.  (In,: 
DO  1/  IP  =  1,N 
CO  12  IC  =  1/N 
l { i / i I C )  =  3. 

jo  1  z  <  -  1 , n 

C(IS,IC)  =  c  (is,  i  c )  +  A  ( i  ?. ,  k  )  *  c  c  s  (  k  ,  i  c: ) 

F  =  F  ♦  1  . 

DU  13  IS  =  1,3 
DO  13  IC  =  1  ,  N 

IF  (AES  (CCI.  (IH,1C)  )  .  IT.  IF  ST*  ASS  (PHI  (I  2, 1C)  )  ) 

3D  IC  50 
C  0  3  I  i  N  "J  £ 

CALL  ?  F  C  E  {IFF  Li,  N,  3,  u, CEL,  NO,  3D, ID) 

CALL  PICE  (TEIL,  B,N  ,N,  3,C2,  SD,i!1D,lD) 

DU  1+  ID  =  1,3 
Cu  14  IC  =  1  ,  '1 

Cl  v  I  i  ,  Ic  )  —  D  lL  (  ~  C  2  (IC,  i.1) 

•  Z  I  J  ci  2 
END 

i  'J  j  DC  r!T  INF  COmiES  »  HE  cmriil  IAIN  “.ALFIX  AN 
I  A  ACE 


J  3  J  :•  U  ■]  1 1  M  C  A  I  3  ( t  .<  K  ,  I  7  Y  3.  1  ,  £  ,  r-  ,  ?  H  I ,  H  ,  N  ,  3 , 0 , 3  I  , 


V.  3  I T  c;  (  8  ,  J-5  JO)  (C11  (F.)  ,K=1,  10) 

-.kite (8,  H'j?.)  (H'iiii;;)  ,.x=i,  U) 
9  fc  j  o  F  ti  a  1  A 1  { L  1  J  .  4 ) 

RITE  (8 , 4  1  C) 

4  10  FOP 91, Vi  ( '  1  •  ) 

4  SITE  (4, 7717) 

A  L  L  ?  L  C  x  P  (  x  I  .i  ;  ,  A  1  1  ,  4  !-i  ,  1 ) 

CALL  e  L  0  I  r-  { 'i  x  9,  n  ,  (i  3  1  ,  L  v1,  ,  -  ) 

4  FITS  (8,4  K) 

L  A 1 1 E  (4  ,  7  7  77) 

CALL  i?LOrr*  (IlJi,  k'H,MAi,  1) 

CALL  ?  L  3 1  ?  Ci:iZ,YHF,LN,3) 

L  L  =  4  n  +  2 

LCN3  (LL-1)  -ICO 
Lc4G(LL)=  187 
LA  *  (LL-1)  =0 
LAT  (LL) =  t 0 
(.BIT  Z  (4, 77  77) 

CALL  PI.CiT  (LONG, LAI, LL,  1) 

CALL  PLCII  (Xu, X H , N N , 2 ) 

CALL  3 1. C  : F  (LC4  j  ,  L  A'l ,  LL  ,  1) 

CALL  ?  L  C  I  ?  (  a  ii  ,  X ,  .'t  'J ,  2  ) 


*****  j p I x £  7  A  1 1 


.  4  i.  C  ,  ij  C  7 


4  UTICZ  c  F  V  A  LUES 

c  it  -  ( i J ,  4 2 c )  m:i,u  ( i ,  1)  ,a  (1, 1) 

fi  FITE  (  4,24  2  a)  44 , 4  ( 1 ,  1 )  ,  F  (  1  ,  1 ) 
2  U  2  J  r  C? .1  A T  (14,  j  X  ,  1  ,V  ( 1  ,  1 )  =  •  ,  F  3.  3 , 3  X  , 


•n  L  Z  £  (4,  77  /  7) 

^172(6,  ^9  2)  ill"!  (I)  ,1  Ari  (I)  ,LOE.G  (I)  ,In  <I)  ,2CE  (I)  ,1  =  1 
1  ,M) 

iVPITE  (o,  39  3)  (El ut  (I)  ,  ELAT  (I)  ,  ELC«G  ( I)  ,  YKMri  (I) 

1  ,  XK.Ito  (I)  ,1=1,  NZ) 

39  i  PCI-MAT  ( 1  X  ,  b  E  6-2) 

« .t  I  -  «  (9  /  /  7  7  7 ) 

*3112(4, 3S3)  ( I  I.'l  3  ( I )  ,  L  A  I  ( I )  , LG  J9  ( I )  ,?.<(!)  ,  X  :i  ( I ) 

1  ,  I  =  1  ,  i!  Z  ) 

«x:r£  (4,  3  9  j)  (BIT  12  (I)  ,  6  LA  I  (x)  ,  jL^iJ  (1)  ,  it  .If  (I) 

1  ,X<.1c  (1)  ,1=1, N7) 

INTI  '!  Z  =  I  M  T  ( 1 1  '  E  ( X )  ) 
p  L  T 1 .1 F  =  B  LCAT  (I  5  11.12) 

Ir  { ?  L 1 1 .1  E  .  5  E  .  x  I  ’1 2  (*)  )  JUIC  j-J«? 

n  l  L  —  —  —  (  —  1.  t  i  —  ,  — ■  4.  v.*  .  *  .3  ,  —  a.  .41,2-  —  1..,.*.  —  ^  ,  —  v,  1  j  ) 

)  9  C  C  N L  I  N  J  2 

***  ..’KITE  ECL 1  A;.D  it  j2  1C  FILE  'EC.  2 .» A  ' 

*:*_.L£(C,J<i2) 

i  9  EG -i.1  A I  (//,  bX,  •  JULIAN  EE  EG  2  1 

1  2032’) 

ITc  ('4,-32  i )  ( J  J  -  r  7  (IEJc)  ,  3921  (I  2Gb)  ,  2<:2Z  (15-LL) 
i  ,  i  3c:j=  i ,  xz) 

Vi  Air  E  (  16,  J26)  (  K  JJL  :i-.  (I  EC  ■:)  ,  F,  G  6  1  ( I  El-)  ,  >.C32  (I".C=) 

1  ,  I  2  C  B  =  1  ,  N  Z ) 

l  :3  5  C  2 .1  A I  ( 1 1  6  ,  z  F  1  5 . 2 ) 

DG  J^7  I=1,JN 

E  i  ( 1 )  =  Y  K  .1  6  ( I )  -  2  L  A  1  ( I ) 

iX  (I)  =  X  ^  .1  a  (I)  -JLONG  (I) 

I  I  IE  ( I)  =  EEC  A  -  (  I ) 


El 

-• 

j  J  •: 

:  H  T 

A  N  i  £ 

*c> 

" •:  1  (1,2)/(FXX1 1  (1, 

1) 

'  ;  • 

32 

.<= 

O  • 

;i\  x 

1(1, 

1) 

♦  FKICI  1  (2,2)  )/4.  +  :0 

- 1 

(1 

/  5 

IN 

*  T 

KE1) 

o  <_ 

Y  = 

(?; 

C  K  X 

1(1, 

D 

+  Pa  X /;  1  (2,  2)  )  /2.  -  PEN 

.1  1 

(1 

IN 

(2. 

,  *  I 

!i  E 1 ) 

sx=  ( (sir,i'<)  **m  5)  *3 

SY=  (  (3132  Y)  **.  5)  * b 
FI=J. 141552c 5/ 12 
CT=CC3  (  T  H  2  1 ) 

3  1  =  31  :■!  ( I  a  1 ) 

DC  19c  1  121119=1,23 
21  =  121  LI F 

X?  (I21LI?)  =SX*COS  (?:*VI)  *CT-SY*J1J  (  =  *; 

IP  (I2LIIF)  =  3  X  *  C  0  3  (  ?  I  *  X I )  *3T*SY*SI.<  *C 

1  4  ;<  1 :  2  ( 1 3  ,  1  932)  Xr-CHLIIF)  ,  ¥  ?  (121.  II?) 

2  ?  3  M  .1 A 1  ( 2  F  1  2 .  7 ) 

*  **  2:0  Or  2  I  1 1  2 3 2  CAlCiil  *\TLC‘-i 


2  2  CONTINUE 


a  j  ;  i  ;  _  X  i',  o  4  >  J  E  c 

l'.?I  22  (5,22  2  1) 

-  -Oi  VO  (//,  •  1*0.72  ( J  1)  ,  ::Oj(JJ),  iK.0  4(.JJ) 
1  ,  IK.Ic(JJ)  ') 

IFnE  (b,  2223)  (TX.-172  (JJ)  ,  1X343  (OJ)  ,  2  K  *1  2  4  (J.J) 

i  ,  j  j  =  1  ) 

«?.  112(3,2225) 

5  AM  (//,  »  IK  P  72  (31),  TXT4-J  (r!v)  ,  '1X224  (1 M  ) 

1  T  .<  2  r>  ( 2  1 )  ’  ) 

-ri:i(b,  222J)  (  F  K  J  7  (31)  ,  1^4c  (3  1)  ,  1XP2  4  (‘O.) 


[•TO 
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